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ABSTRACT 

Using proper motion data for 894 stars in the Orion Nebula Cluster (ONC) compiled 
by Jones & Walker in 1988, we search for binaries with apparent separations in the 
range 1000-5000 AU, and find an upper limit of three. Using a Monte Carlo method, 
we test the consistency of this result with two hypotheses: i) that the cluster contains a 
binary population identical to that found in the solar neighbourhood, and ii) that the 
cluster contains no binaries at all in this separation range. We obtain results strongly 
favouring the latter hypothesis. 

Star formation in the Galaxy is seen to occur in a variety of different environments, 
but it has been proposed that most stars may be formed in dense regions similar to 
the ONC, rather than in less dense groupings like that found in Taurus- Auriga. Since 
roughly 15 per cent of galactic field stars are known to be in binaries with separations 
greater than 1000 AU, the apparent absence of such binaries in the ONC places an 
upper limit on the contribution that dense clusters can make to galactic star formation. 

Key words: stars: formation - stars: pre-main-sequence - stars: statistics - open 
clusters and associations: individual: Orion Nebula Cluster - binaries: general. 



1 INTRODUCTION 



198C ), contains so me 1000 stars with a peak densi t y of 2- 
5xl0 4 stars pc~ 3 (|McCaughrean fc Stauffer 1994 



Hillen- 



Among st the stars around us in the Galactic Held, single brand fc Hartmann 1998| ). This compares to core densities 



systems like our own are a minority. Most stars are tound 
in binaries, or in systems of even higher mult iplicity. Re- 
cent surveys o f the nearby stellar pop ulation (Duquennoy 



|fc Mayor 199l| [Fischer fc Marcy 1992 ), extending to 20 pc 



or so, indicate that about 60 per cent of systems in the so- 
lar neighbourhood are binary. In younger populations, such 
as those in nearby star-forming regions, the degree of bi- 
narity is often even higher. Surveys of the low-mass dark 
cloud complexes in Taurus-Auriga, Ophiuchus, Chameleon, 
Lupus, and Corona Australis reveal a fraction of binary sys- 
tems roughly twice as high as that of the solar neighbour- 
hood in the ra nge of binary sepa r ations detected (typicall y 



10-1500 AU) flGhez et al. 1997| ; |Kohler fc Leinert 1998| ). 
Thus single stars are rare in these clusters. An exception to 
this trend is the Orion Nebula Cluster (ONC) in the Orion A 
giant molecular cloud, which has a binary fraction roughly 
the same as the solar neighbourho od value, at least over the 
separation range of ~ 25-500 AU (Prosser et al. 1994 ; 



Close & Beck 1999) 



gctt, Strom fc Chez 199^ ; |Petr et al. 199$ |Petr 1998} 



Pad- 



in the Taurus- Auriga complex of only around 10 stars pc~ 3 . 
It has been argued, based on estimates of star formation 
rates ( Miller fc Scalo 1978 ), that dense OB associations like 
the ONC are responsible for the majority of star formation 
in the Galaxy, and are therefore representative of the typical 
environment for early stellar evolution. 

The results on binarity in star-forming regions would 
seem to support this conclusion. However, the high density 
of the ONC has meant that the surveys carried out there 
so far have only been able to detect binaries separated by 
less than about 500 AU. It would be very useful to examine 
the binary fraction at wider separations, but in doing so one 
is increasingly likely to observe chance projected alignments 
rather than actual bound pairs. Bate, Clarke & McCaugh- 
rean (1998) analysed the mean surface density of compan- 
ionsin the centre of the ONC us ing da ta from Prosser et al. 
(1994) and McCaughrean et al. (1996), and found weak evi- 



The ONC differs from the other regions mentioned in 
that it is densely clustered, and contains a number of mas- 
sive O and B type stars. In particular, the core of the ONC, 
known as the Trapezium cluster after the eponymous four 
central OB stars which dominate it (Herbig fc Terndrup 



dence there for a deficit of binaries with separations greater 
than 500 AU. But to do better, and to detect wider binaries 
properly, one must use the fact that they have a common 
proper motion across the sky. Specifically, one looks for stars 
whose relative velocity is less than the critical value needed 
to escape their mutual gravitational attraction (which can 
be determined by making an assumption or estimate of the 
stars' masses). 
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In this work we attempt to detect common proper mo- 
tion binaries in the ONC using the data obtained and cat- 
alogued by Jones & Walker (1988; hereafter also JW), and 
thence determine whether the number found is consistent 
with the field star binary distribution over the same separa- 
tion range. 



2 FINDING COMMON PROPER MOTION 
PAIRS 

The JW study used a series of optical wavelength photo- 
graphic plates taken over a baseline of about 20 years, and 
obtained proper motion measurements for over 1000 stars. 
Compared to deeper and higher spatial resolution surveys 
carried out using the Hubble Sp ace Telescope at optical 
wavelengths (Prosser et al. 1994) and ground-b ased adap- 



tive optics techniques at infrared wavelengths (VlcCaugh 



rean & Stauffer 1994), the JW data reveal only a fraction 



of the total stellar population. However, it remains the best 
proper motion survey to date, and covers a relatively large 
area on the sky, extending roughly 15 arcmin or 2 pc from 
the cluster centre, taken to be 6 1 Ori C, the most massive 
of the four Trapezium stars. 

Ideally, one would have data in three dimensions for po- 
sition and velocity, plus the mass of each star, as it would 
then be a simple task to identify all the bound pairs within 
the cluster, and give an accurate value for the binary frac- 
tion. Instead, because the data provides information only in 
two dimensions, one can do no more than place an upper 
limit on the number of possible bound stars. The approach 
taken here is to look for apparent binaries in the JW data 
- pairs whose 2D positions and velocities allow them to be 
bound if one ignores the third dimension - and then to com- 
pare the number of these with the number of similar pairs 
observed in a randomly generated model cluster with known 
parameters. 



Hillenbrand (1997) found the mean stellar mass in the 
Trapezium to be approximately 0.8 Mq; here we make the 
assumption that all the stars are of solar mass. The distance 
to the ONC is taken to be 470 pc (the value used by JW), 
at which 1 pc corresponds to 7.3 arcmin on the sky. We test 
for binarity at projected separations up to 5000 AU, beyond 
which the critical relative velocity (the escape velocity for a 
pair at that separation) is lower than the mean error in the 
JW velocity data (~ 0.8 kms -1 ). A lower limit of 1000 AU is 
imposed by the observational resolution of ~ 2 arcsec. Thus 
our study is complementary to previous work, where binaries 
were sought at separations roughly ten times smaller. 

An analysis of the JW data gives the following initial 
results. Eliminating stars with less than 80 per cent mem- 
bership probability leaves 894 stars in the cluster catalogue, 
which give a total of 894 x (894 - l)/2 = 399171 pairs. Of 
these, 192 have apparent separations between 1000 and 5000 
AU. Assuming solar masses and using the JW velocity mea- 
surements (ignoring errors for the time being) , three of these 
192 have relative velocities lower than the critical velocity 
for their separation. 

It is perhaps likely that if data for the radial dimension 
were known, these binaries would turn out to be unbound. 
But it is also likely that the effect of errors has been to 
'disrupt' other apparent binaries in the cluster. These com- 
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Figure 1. The radial density distribution of a typical model clus- 
ter (x symbols) comp ared to that of the ONC as surveyed by 
Jones & Walker ( 198a ) (triangles). The solid curve shows what 
gradient would be expected for a singular isothermal sphere of 
radius 5 pc. The effect of the finite boundary can be seen as a 
slight downturn with increasing radius. Similarly, the effect of the 
flat core can be seen in the data points as a flattening towards 
small r. The dashed line marks the core radius. 



plications - in particular the projection effect - mean that 
we cannot directly compare our result with a prediction from 
the Duquennoy & Mayor (or any other) period distribution. 
Instead, as mentioned above, the comparison is made by 
constructing an ensemble of model clusters whose observ- 
able parameters match the JW data statistically, but which 
contain the population of binaries we want to compare with. 



3 SIMULATED CLUSTERS 

The model clusters have a three-dimensional density distri- 
bution corresponding to an isothermal sphere with a flat 
core: 



p(r) 



P0 



if r < R c 
iir>R c 



(1) 



where R CO rc = 0.04 pc. A similar distr ibution was shown by 
Bate, Clarke, & McCaughrean (199S ) to give a good fit to 
the projected density distribution of the JW data. Figure 
fy shows a comparison between the simulated and observed 
data. A finite radius of 5 pc is chosen, within which a suf- 
ficient number of stars are distributed according to (hi) to 
ensure that, after the corrections discussed below, roughly 
900 appear projected within the same radius as the JW sur- 
vey (~ 2 pc). Each system is then assigned a velocity chosen 
from a Maxwellian distribution with dispersion of 4 kms , 
equal (in three dimensions) to that found by Jones & Walker 
for most of the cluster. 

Two different binary populations are considered, one 
with a binary fraction of zero (i.e. no binaries at all), and the 
other with a population matching that found by Duquennoy 
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& Mayor (1991). In the latter, the periods have a log-normal 
distribution peaking at 180 years, which corresponds to a 
typical separation for solar mass stars of 30 AU. 

Once a cluster has been generated with these parame- 
ters, we mimic the process of observing it, applying the same 
errors and restrictions to the data as were present in Jones & 
Walker's survey. For example, the velocity errors in the JW 
catalogue were found to be log-normally distributed with a 
peak at 0.8 kms -1 , and we apply noise to the velocities in 
the model cluster using the same distribution. Similarly, the 
JW survey had a resolution limit of about 2 arcsec or 1000 
AU, so any model stars which appear closer in projection 
are merged into one. 

We must also correct for the fact that some of the bi- 
naries present in Duquennoy & Mayor's population would 
have appeared as single stars to Jones & Walker, because 
one companion would have been fainter than they could de- 
tect. In the outer part of the cluster, JW give their detec- 
tion limit as / ~ 16, but in the centre, where there is bright 
background emmission from the Orion Nebula, this limit 
is reduced. Indeed, within about 2 arcmin (0.3 pc) of the 
centre they apparently detect no stars fainter than I = 14. 
We take the mean reddening of Ay — 2.4 magnitudes mea- 
sured towar ds the optically-visible stars by Herbig & Tern- 
drup (1986), which corresponds to a mean reddening in the 
I band of 1.4 magnitudes. Then, in the outer regions, the 
JW study is sensitive to stars with intrinsic / magnitudes 
as faint as 14.6, which at the assumed distance and an age 
of 1 Myr (see §4) corresponds to a mass of roughly 0.2 Mq 
(b'Antona & Mazzitelli 19941) . The brighter detection limit 



within 0.3 pc corresponds to a mass of 0.6 Mq. 

The minimum companion mass in the wide binaries 
included in Duquennoy & Mayor's survey was also about 
0.2 Mq, this limit being determined by the magni tude range 
of the data they used - see Patience et al. ( 1998 ) . We ther- 



fore choose to apply a correction only to the centre of our 
simulated clusters, removing a certain fraction of the binary 
components whose projected distance from 9 1 Ori C is less 
than 0.3 pc. Denoting this fraction by X, we need to evaluate 
the following integral: 



X 



Jp.2 s ( m ,^) d ^ 



f(m) dm 



(2) 



where f(m) is the stellar mass function in the cluster, and 
s(m, fj,) is the companion mass function - the distribution of 
companions of mass /j, to primaries of mass m. For simplicity, 
we assume that the two mass functions match in the region 
where //<///: 



s(m,fi) 



/(/i) if /i < m 
if /i > m 



(3) 



(The pri mary i s defined to be the more massive star.) Hil- 
lenbrand (1997) found the mass distribution of stars in the 



ONC to be reasonably fitted by a Miller-Scalo mass function 
- a semi log-normal, falling from a maximum at around 0.1 
Mq. Using this function for f(m), we evaluate the integral 
numerically, and obtain a result of 0.75 for X. 



With these corrections applied, each model cluster is then 
analysed for apparent binaries, as was done with the JW 
data in the previous section. The results for ensembles of 



400 - 



300 - 



S»200 - 



100 - 






10 15 

Apparent binaries 



20 



25 



Figure 2. Histogram showing the results of observing three en- 
sembles of simulated data, one having no binaries, and the other 
two (hatched) having a binary population matching that found in 
the solar neighbourhood. Each ensemble comprises 2000 clusters. 
The data with heavy hatching has been corrected for a brighter 
detection limit in the cluster centre than in the outer region, 
whereas the data with light hatching assumes the fainter detec- 
tion limit throughout (see text, §3). The dashed line indicates the 
observational result of three apparent binaries for the ONC. 



2000 clusters, plotted in Figure bl show that the clusters 
with zero binary fraction, when seen in projection, are far 
more likely to show three common proper mo tion binaries 
than the clusters with a Duquennoy & Mayor (1991) binary 
population. Specifically, we find that in the former case, a 
result of three apparent binaries occurs in 427 (21 per cent) 
of the 2000 clusters, while for the Duquennoy & Mayor clus- 
ters which have been corrected for mass detection limits it 
occurs in just 13 (0.7 per cent) of them. A result of three 
binaries is 49 centiles away from the median of the latter 
distribution, which in a gaussian curve would correspond to 
a distance of just over 2.3 times the standard deviation. 



4 DISCUSSION 

Over the bulk of the cluster, Jones & Walker found a 
roughly constant three-dimensional velocity dispersion of 
about 4 kms' 1 , increasing only slightly towards the core. 
They also found some evidence for anisotropy in the velocity 
dispersion towards the outer regions of the cluster, increas- 
ing in the radial direction - the direction towards and away 
from the cluster centre - and decreasing in the tangential 
direction. In terms of separation, 4 kms -1 corresponds to 
a hard/soft binary limit of about 30 AU, much closer than 
the binaries we have considered here. Our wide binaries are 
therefore well into the regime where one would expect dis- 
ruption through encounters with other cluster members. 

The typical encounter time r of a binary with another 
star in the cluster can be crudely estimated as 
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l 

naS 



(4) 



where n is the stellar density, a the velocity dispersion, and S 
the binary cross section (taken simply to be ira 2 , with a the 
semi-major axis, since we ignore the effect of gravitational 
focusing for wide binaries). In the core, taking a stellar den- 
sity of 2 x 10 4 pc -3 , we expect a binary with a = 1000 AU 
to have r less than 10 years. Further out, r increases as the 
density drops, surpassing 10 7 years at about 1 pc from the 
cluster centre (assuming an isothermal sphere distribution). 
In her s tudy of masses and ages in the ONC, Hillen- 
brand ( 1997 ) found the great majority of stars to be less 
than a million years old. Since this is significantly shorter 
than the disruption timescale for wide (1000-5000 AU) bina- 
ries in the outer regions of the cluster, their apparent absence 
could be taken either as evidence that none were formed in 
that separation range, or as an indication that all the stars in 
the outer cluster have at some stage passed through a denser 
environment, such as the core. Taking the core radiu s to b e 
about 0.2 pcjj as found by Hillenbrand & Hartmann (1998), 



a star moving at 4 kms~ will take some 10 years to pass 
through, so that even one visit would suffice to disrupt most 
of these wide binaries. But the ONC is a young cluster, and 
it is not clear that the majority of stars will have visited the 
core in the time available. 
Brandner & Kohler 



199£) have examined the period 



distribution of binaries in two subgroups of the Scorpius- 
Centaurus OB association, and found significantly more 
wide (a > 200 AU) binaries in the subgroup containing 
fewer massive B type stars. It is plausible that the binary 
period distribution in star-forming regions might be deter- 
mined at the time of formation, and be influenced by condi- 
tions like the temperature and densit y of th e cloud, as has 
been proposed by Durisen & Sterzik (1994). However, it is 



difficult to exclude the effects of subsequent dynamical evo- 
lution, which can happen on a short timescale in these dense 
environments, and which can account for a deficit of wide 
binaries independently of the initial distribution. This ef- 
fect has recently been investigated in dynamical simulations 
of clusters in equilibrium, cold collap se, and expansion by 
Kroupa, Petr & McCaughrean (199£). An overview of dy- 
namical pr ocesse s in young clusters is provided by Bonnell 
& Kroupa ( |l999| ). 

In conclusion, we observe that about twenty per cent of 
the binaries in Duquennoy & Mayor's survey are wider than 
1000 AU, so with a total binary fraction of 0.6, they com- 
prise around fifteen per cent of G dwarf stars in the Galactic 
field. A similar figur e can be assumed for the M dwarfs in 
Fischer & Marcy's (1992) survey. If such systems do not 
come from regions like Orion - whether because they are 
never formed or because they are dynamically disrupted - 
then a significant number of field stars must be formed else- 
where, possibly in cooler or less dense environments. This 
raises the more general question of what 'mix' of star form- 
ing regions is required to produce the field distribution we 
see today, and demonstrates the utility of binary stars as 



a diagnostic probe for such syntheses. Further collation of 
binary statistics in different star forming regions, over the 
widest range of separations, should continue to shed light on 
these questions. 



5 ACKNOWLEDGEMENTS 

We thank Ian Bonnell, Pavel Kroupa, and Melvyn Davies 
for valuable help and discussions. A. Scally is grateful for 
the support of a European Union Marie Curie Fellowship. 



REFERENCES 

Bate M. R., Clarke C. J., McCaughrean M. J., 1998, MNRAS, 

297, 1163 
Bonnell I., Kroupa P., 1999, in McCaughrean M. J., Burkert A., 

eds, The Orion Complex Revisited, ASP Conference Series, in 

preparation. 
Brandner W., Kohler R., 1998, ApJ, 499, L79 
D'Antona F., Mazzitelli I., 1994, ApJS 90, 467 
Duquennoy A., Mayor M., 1991, A&A 248, 485 
Durisen R. H., Sterzik M. F., 1994, A&A 286, 84 
Fischer D. A., Marcy G. W., 1992, ApJ, 396, 178 
Chez A. M., McCarthy D. W., Patience J., Beck T., 1997, ApJ, 

481, 378 
Herbig G. H. & Terndrup D. M., 1986, ApJ, 307, 609 
Hillenbrand L. A., 1997, AJ, 113, 1733 

Hillenbrand L. A. & Hartmann L. W., 1998, ApJ, 492, 540 
Jones B. F., Walker M. F., 1988, AJ, 95, 1755 
Kohler R., Leinert Ch., 1998, A&A, 331, 977 
Kroupa P., Petr M. G., McCaughrean M. J., 1999, MNRAS, in 

preparation. 
McCaughrean M. J., Rayner J., Zinnecker H, Stauffer J., 1996, 

in Beckwith S. et al., eds, Disks and Outflows Around Young 

Stars, Lecture Notes in Physics 465. Springer, Heidelberg, 

p. 33 
McCaughrean M. J., Stauffer J. R., 1994, AJ, 104, 1382 
Miller G. E., Scalo J. M., 1978, PASP, 90, 506 
Padgett D. L., Strom S. E., Ghez A. M., 1997, ApJ, 477, 705 
Patience J., Ghez A. M., Reid I. N., Weinberger A. J., Matthews 

K, 1998, AJ, 115, 1972 
Petr M. G., 1998, Ph.D. thesis, University of Heidelberg 
Petr M. G., du Foresto V. C, Beckwith S. V. W., Richichi A., 

McCaughrean M. J., 1998, ApJ, 500, 825 
Prosser C. F., Stauffer J. R., Hartmann L., Soderblom D. R., 

Jones B. F., Werner M. W., McCaughrean M. J., 1994, ApJ, 

421, 517 
Simon M., Close L. M., Beck T. L., 1999, ApJ, in press. 



Note that this is not the same core radius as was used pre- 
viously in generating the model clusters. That value (0.04 pc) 
simply parameterised the isothermal sphere distribution used to 
match Jones & Walker's data. 
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